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In this work we investigate the contribution to Bd, s mixing from both scalar and vector unparticles 
in a number of scenarios. The emphasis of this work is to show the impact of the recently discovered 
3a evidence for new physics found in the CP phase of B s mixing. Here we show that the inclusion 
of the CP phase constraints for both Bd and B s mixing improves the bounds set on the unparticle 
couplings by a factor of 2 ~ 4, and one particular scenario of scalar unparticles is found to be 
excluded by the 3<r measurement of <j> s . 
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CN I. INTRODUCTION 

Cd It has recently been suggested [TJ |5] that there may exist a nontrivial scale invariant sector at high energies, known 
* as unparticle stuff. These new fields with an infrared fixed point are called Banks-Zaks fields [3], interacting with 
O^l standard model fields via heavy particle exchange, 

Oh 

■ Here Osm is a standard model (SM) operator of mass dimension dsM, Obz is a Banks- Zaks(BZ) operator of mass 
dimension dgz, with k = d$M + dsz — 4, and My is the mass of the heavy particles mediating the interaction. At a 
(-H scale denoted by Au the BZ operators match onto unparticle operators with a new set of interactions, 

\d-BZ-du 

m C u " SM Ou (2) 

> M u 
O 

l^-} where Ou is an unparticle operator with scaling dimension du and Cu is the coefficient of the low-energy theory. 
£f~) Unparticle stuff of scaling dimension du looks like a nonintegral number du of invisible massless particles. In this 
work we shall study both scalar unparticles {Ou) and vector unparticles {O u ) with couplings to the SM quarks as 
follows, 
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S,q'q S,q'q 



Scalar unparticles: .^^(i _ <^ + ja ^(i + ^fl/* ^ (3) 



y< q'q „ v > q'q 



^ Vector unparticles: ^— ^ q'j^l - 75 )g O u + q'j^l + j 5 )q O u (4) 
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Here we assume that the left-handed and right-handed flavour-dependent dimensionless couplings c£, c s R , and 
are independent parameters. We analyze a number of scenarios, in each case determining the allowed parameter space 
and placing bounds on the unparticle couplings. 

The propagators for scalar and vector unparticle fields are as follows j2j H] , 

d*xe iP '*{0\TO u (x)O u (m = * 9 Ad " (p2 J^ du ^ u (5) 

2 sin dyTT (P z 4- te) z au 

d 4 xe lP -(O|T^(x)0^(O)|O) = i . d " 1 9 (p2 \.y 2 _i (6) 

2sina^7r \P* + ie) z au 
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where 

16^ 5 / 2 T(d u + 1/2) 

Since the publication of the first theoretical papers on this new subject [TJ [5] there has been a huge interest in 
unparticle phenomenology, for example, Collider signatures |ll3[i[ail[S[ia[ni[Il[I3[Tl[I3[Tni[ni[Tl[ini 

EH1 EH E21 E31 EH EH Ell E71 EH Ell Eol EH E21 EH EH El], CP violation [35J E3 EH EH Sol E], meson mixing 

[H ESI HH H31 HH HS1 HE1 HZI. lepton flavour violation gH S3 EH1 EH [H |S1], consequences in astrophysics 
[551 1551 1571 1551 1551 1501 I5T1 [55] , in neutrino physics [531 IM1 [55J [55J [57] and in supersymmetry [551 ESI [70] • In this work 
we shall study the constraints coming from the measurements of B s d meson mass differences AM S ^ and also their 
CP violating phases (f> s .d- In the Bd system these quantities have been well measured for some time and show only 
small deviations from the SM expectation. In the B s system recent measurements have also found small discrepancies 
between the SM expectation for AM S [71], but now the CP violating phase cj) s , measured by the D0j75] and CDF[75] 
collaborations reveals a deviation of 3<r[7U [75J [53]. This is the first evidence for new physics in b «-> s transitions. 
Studying both vector and scalar unparticles, we study the constraints imposed by these latest measurements on the 
coupling between SM fields and unparticles, with particular interest on the impact of <f> 3 . 

In Sec. ITT] we discuss unparticle contributions to general meson-antimeson mixing and in particular, to the case of 



B mixing. Section III contains the results of our numerical analysis where we first consider scalar and then vector 



unparticle effects in Bd s mixing. In each case we set bounds on the allowed parameter space using constraints from 



measurements of AM djS and from <p d _ s . Finally, Sec. IV concludes our results 



II. MESON-ANTIMESON MIXING FROM UNPARTICLES 



With unparticle operators coupling to standard sodel operators as in Eq. 
to contribute to meson-antimeson mixing via the s- and i-channel processes shown in Fig. [1 



3|4) it is possible for unparticle physics 




FIG. 1: s- and t-channel unparticle contributions to meson mixing. 

Using the interactions listed in Eq. ( 3p l to evaluate the s- and t-channel contributions to meson mixing as shown 
in Fig. [JJwe obtain the effective Hamiltonian for scalar unparticles as, 



H 
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eff 
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m q c R ~ m q'C L 



S.q'q 
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! Qi (m q >c S £ q q - m q c^ q q ^j ( 
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m q c L q - m q , C ji q 
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and for vector unparticles we obtain, 



V,q'q _ 
off — 



A 



du 



2sind W 7r Aff u ~ 2 
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f3-d u s 3-d u 



V,q'q V,q'q 

m q c H — m a > c 



1 



tf-du 



c2 — du 



Q 2 (m q > cl' q q - m q c^ q 9 ) ( 

i A / V.q'q V.q' q\ ( V.q 'q V.q q\ 

+ Q2 [m ql c£ - m q c L ' H J \m q c L H - m q < c H '* H \ 
+ 2 Q 4 (m q > c^ q q - m q 

'1 (c^y+Q, (c v / 



V.q'q V.q'q 

m q c L ' q - m q , Cji q 



(9) 



4Q 5 (c V / 



q „S,q'q 
C R 



In the final line we have used a Fierz identity to rearrange the operator (V — A) ® (V + A) into the scalar operator 
Q5. Here H.^ s describes the effective Hamiltonian for the case of scalar unparticles and TC^ S for vector unparticles, 
hence the total effective Hamiltonian is simply, 



T-ll'l — n-/S,q'q 
n cff — n efi 



V,q'q 
eff 



Above we have defined the quark operators Q\ — Q5 as follows, 



Writing the effective Hamiltonian in terms of these operators we have, 

n qq = £(Cf + CY)Q l + £(0? + Cl)Qi 



(10) 
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= qTwZ <tM L 


(11) 
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= q'rQI Jfgfi 


(12) 


Qs 


—/a 8 —10 a 

= QrQl QrQl 
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Qi 


—l a a —iP 8 

= Qr Ql Ql Qr 


(14) 


Qs 


—/a fj _//3 a 

= Qr Ql Ql Qr 


(15) 



(16) 



Here the operators Qi,2,3 are obtained from Qi,2,3 by the exchange L <-> R. 

The hadronic matrix elements, taking into account for renormalization effects, are defined as follows, 

1 



(M°|Q 2 ( At )|M°) = 

(m°|Q3(m)|m°> = 

(m°\q 4 ^)\m ) = 

(M°\Q 5 (H)\M°) = 

where, i?(/i) 



M M f 2 M R^)B 2 ( f i) 



24 

^M M f 2 M R^)B 3 ^) 
\M M f M R^)BM 



M M f 2 M R^)B 5 ^) 



(17) 
(18) 
(19) 
(20) 
(21) 
(22) 



M M 
m q (n) + m 9 /( / u) 

From Eq. ([&||9| , it is straightforward to calculate the Wilson coefficients for scalar unparticles which are as follows 

du , 



c* 



A 



e -i<f>u f M M 



smduK Mi r \ Af, 



< S.q'q S.q 

m q ,c L qq -m q c R q 
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du 



e -i<Pu fM M 



C> = 



sin dun Mf 
A d, 
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S,q'q 
m q' C L m 



I dun Mt 



e -i<ht ( M 2 M 



M 



du 



S.q'q S.q 
™q' C R ~ m q C L 



^ ( m q C R q q 
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S.q'q S.q'q 
m q C L ~ m Q' C R 



S.q'q S.q 

m q c L q - m q >c^ q 
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°5 — °1 — °3 



(23) 

(24) 

(25) 
(26) 
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and for vector unparticles we find the following Wilson coefficients, 




CY = = (33) 

where we have approximated t — s ~ M M . From these two sets of Wilson coefficients it is clear that the case of 
vector unparticles not only includes more contributions as C\ ^ 0, C\ ^ 0, but also that their Wilson coefficients 
are enhanced by a factor Ku/M M compared to the scalar unparticle case. As a result the vector unparticle parameter 
space shall be suppressed by the same factor. 

These Wilson coefficients will mix with each other as a result of renormalization group(RG) running down to the 
scale of Mm- For the B system, with a scale of new physics Ay = 1 TeV, these Wilson coefficients at the scale 
Hb = mb may be approximated as [77] . 



Cifa) « 0.805 Ci (A«) (34) 

C 2 (nb) « 1.988 C 2 (Au)- 0.417 C 3 (A W ) (35) 

CM « -0.024C 2 (A W ) + 0.496 C 3 (A U ) (36) 

CM « 3.095 C 4 (A U ) + 0.725 C 5 (A W ) (37) 

C 5 ( Mb ) « 0.086 C 4 (A W )+ 0.884 C 5 (A U ) (38) 



The AF = 2 transitions are defined as, 

(MV C A /= 2 |A/ ) = M 12 (39) 

with the meson mass eigenstate difference defined as, 

AM ee M H - M L = 2\M 12 \ (40) 

We can define in a model independent way the contribution to meson mixings in the presence of new physics (NP) 
as, 

M12 = M! S 2 M (1 + R) (41) 

where Mf 2 M denotes the SM contribution and Ree re" 7 = M^ p /M™ parameterizes the NP contribution. 
The associated CP phase may then be defined as, 

<i> EE arg(M 12 ) = SM + NP (42) 

where SM = arg(M^) and NP = arg(l + re ia ). 

At this point it is important to make a further remark regarding the unparticle like model derived from a scale 
invariant theory of continuous mass fields as discussed in [7S]. Although this model contains no fixed point or 
dimensional transmutation and therefore does not correspond to unparticles as defined in [2], it does contain unparticle- 
likc local operators coupling to the SM. In this model it can be shown that the scale invariance will be broken by 
interactions with SM fields, resulting in a mass gap which could be rather large |79j . In this setting we may then 
expect that there are unparticle states lying below the scale of 1 TeV. Therefore the above approach of integrating out 



5 



unparticle effects at the scale of 1 TeV would not be valid in this special case, rather the Wilson coefficients should be 
instead directly calculated at the low-energy scale. Such a calculation would offer its own challenges, in particular the 
i-channel exchange would involve a nonlocal hadronic matrix element. In this work we shall not consider such a case 
and instead assume that our unparticles do not derive from a scale invariant theory of continuous mass fields. Rather, 
we follow closely along the general framework set out by Georgi's original work in [2] where there is dimensional 
transmutation at the scale A^. 



A. Bd, a mixing and unparticle physics 

In this work we shall focus on the constraints imposed on unparticle physics couplings from Bd. s mixing. Therefore 
we set q' = b, q= s,d and M° = B°, B°. 

In the B system, the standard model contribution to M\ 2 is given by, 

Mli U = ^f M B q V B f Bq B Bq (v; q V tb f S (x t ) (43) 

where Gf is Fermi's constant, M\v the mass of the W boson, fj B — 0.551 80J is a short-distance QCD correction 
identical for both the B s and Bd systems. The bag parameter Bg q and decay constant fg are nonperturbative 
quantities and contain the majority of the theoretical uncertainty. Vt q and Vtb are elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [53152], and S {x t = fhf/M^) = 2.34±0.03, with m t (m t ) = 164.5 ±1.1 GeV [S3], is one of 
the Inami-Lim functions |84j . 

We can now constrain both the magnitude and phase of the NP contribution, r q and a q , through the comparison 



of the experimental measurements with SM expectations. From the definition of Eq. (41) we have the constraint 



AM, 



Pi = AJV-^pf = V 1 + 2r " C ° S °« + r " (44) 
The values for p q given by the UTfit analysis [7H[7B] at the 95% C.L. are, 

Pd = [0.53, 2.05] (45) 
Ps = [0.62, 1.93] (46) 

These constraints on p q encode the CP conserving measurements of AM^ s . The phase associated with NP can also 
be written in terms of r q and a q , 

■ .NP r q smcr q 
sm<p q = ■ : 



1 + 2r q cos a q + r q 



cos0 np = , l + r,CQ6 gg (4?) 
1 + 2r q cos o q + r 2 



Here [ZH [76] gives the 95% C.L. constraints, 

0NP = [_ 16 . 6j 3 .2]° (48) 
0NP = [-156.90, -106.40]° U [-60.9, -18.58]° (49) 

these constraint represent those of the CP phase measurements of 4>d,s- 

In order to consistently apply the above constraints all input parameters are chosen to match those used in the 
analysis of the UTfit group [74j [76] with the nonperturbative parameters, 



fB s \/Bg 3 = 262 ±35MeV (50) 

£ = 1.23 ± 0.06 MeV (51) 

fg 3 = 230 ±30 MeV (52) 

f Bd = 189 ±27 MeV (53) 



G 



III. NUMERICAL ANALYSIS 

In our analysis we shall first consider the case of scalar and vector unparticles separately and further subdivide each 
into four classes as follows, 

• One real coupling; cl ^ and cr = 0, with £ K 

• Two real couplings; cl ^ and cr ^ 0, with {cl, cr} € K 

• One complex coupling; cl 7^ and cr = 0, with cl € C 

• Two complex couplings; cl ^ and cr 7^ 0, with {cl, cr} G C 

Hence we have a total of eight scenarios to consider in the following analysis. 

For our analysis we shall take the unparticle scale A u = 1 TcV. In the literature 0[36ll!2l!3IHimi46l|47] it 
has been common to perform an analysis with a fixed value of the scaling dimension, du = |. It was recently shown 
[551 [5S] that from unitarity considerations [57] there exists a bound on the scaling dimension, du > j\ + j% + 2 — <5j lj2 .o- 
Here are the operators Lorentz spins. In particular this leads to the bounds du > 2(3) for our scalar (vector) 

unparticle operators, which conflicts with the choice §. First it should be mentioned that these bounds are derived 
from a conformal field theory(CFT) point of view. Here we are considering a scale invariant sector at high energies, 
and in general scale invariant theories need not also be invariant under the conformal group. The specific example of 
Banks-Zaks fields though, is in fact invariant under the full conformal group [551 155] , 




1 1 1 i ' 1 1 1 — ' — i — ' — i — i — ' — i — 

-1 -a.5 a 0.5 1 -4 -3 -2 -1 a 1 2 3 4 
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FIG. 2: Constraints on the du versus cf parameter space from Bd mixing(left) and B s mixing(right) for the case of a single 
real coupling cf 7^ and cf^ = 0. Black points indicate the AMd, s allowed regions, while grey points indicate the regions are 
in agreement with both AM, and the CP phase <f> q . 

In order to make a connection with the literature, and to quantify the impact of the recent measurement of the CP 
phase <t> s , we choose to fix the scaling dimension at du = |. Through simple scaling of the coefficient Sd u , it is also 
possible to make a connection to a value of the scaling dimension in the region suggested by the above mentioned 
unitarity bound. In general the bounds on unparticle couplings with the scaling dimension in the region du > 2(3) 
shall be much weaker than for du — |, unless the value is close to an integer, where there is a pole in the unparticle 
propagator. For example, if we write the scaling dimension for scalar unparticles as du = 2 + 65, with eg a small 
positive number, then we have |<S'd M= 3/2| = \Sd u =2+e s I f° r a value of £5 w 0.00021. Hence the bounds derived here from 
the CP conserving quantities AM^s shall also apply to the case du ~ 2.00021. The situation is not so straightforward 
for the bounds derived from CP violating constraints of (f>d,s- For vector unparticles we have a similar situation for 
du = 3 + ey, with ey ~ 2 x 10~ n . The above unitarity bound shall be considered in more detail elsewhere [55] . 

A. Scalar unparticles 

In this section we shall analyze the contribution to Bd, s mixing from scalar unparticles taking a number of example 
scenarios. In each case we shall use the measurements of AM^ S as well as the CP phases <fid,s to constrain the 
unparticle parameter space. 
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One real coupling: cf 7^ 0, = 



In this case the Wilson coefficients at the sale Ay simplify to 



c! 



smdifn M M 



A 2 



du 



A 



du 



(M 2 



M 



J sin d u ir M M \ K u 



du 



m m 



(4 M ) 

( S,i 



sind^7r M M 



Ml 
A 2 



CS = Ci 



du 







SM 



(54) 

(55) 

(56) 
(57) 



where q = d, s correspond to Bd or B s mixing. 

In this first case, in order to have two free parameters, we shall allow the scaling dimension du to vary. In general 
the unparticlc effects shall be smaller for larger du- In this case the scaling dimension du also uniquely determines 
the NP CP phase and as such the small CP phase allowed by (f)d and the large CP phase required by the recent 
measurement of <j> s will act as a stringent constraint on the du ~ cf parameter space. 

Fig. ^displays the constraint on the du — cf parameter space from the Bd mixing(left panel) and B s mixing(right 
panel), including CP phase constraints. As expected the general feature of these plots is that for small du the mixing 
contribution from unparticles is large and so the couplings c L ' bq are strongly constrained. For larger values of du the 
opposite is true with du > 2 resulting in no bound on c S ^ bq from either AMd or AM S . 



du 




-* i.* £ £ io- N io- n 10-'" 10- 

" 4 2 4 



FIG. 3: Plot of the variation of the phase and magnitude of the quantity Sd u with the scaling dimension du- 



In the case of scalar unparticles the Wilson coefficients have the common quantity, 



S,i 



A 



du 



sin dun Mm 



du 



(58) 



Here the NP CP phase is determined solely by &vg(Sd u ) such that a q = aig(S du /Mf^ A ). The magnitude of Sd u also 
determines the size of the unparticle contribution to B mixing and hence determines the level of constraint on the 
coupling cf . Fig. [3] shows the variation with du of the phase and magnitude of Sd u , from which we can see that the 
phase of Sd u is constrained to be in the range (— n, 0). The magnitude of Sd u is generally decreasing for increasing 
du except for when du approaches 2 where there is a singularity due to the factor l/sind^7r. The general features of 
Fig. [2] are now easy to understand. The large value of Sd u for small du restricts the allowed values of cf to a very 



narrow region. For example, for du = I we have the bounds 



|cf' M | < 0.067 

I S,bs\ 

\c L I 



< 0.23 



(59) 
(60) 
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The difference in the level of the above constraints on the couplings c L bd and c L ' bs are due to the ratio of CKM matrix 
elements Vtd/Vt s ~ A. Increasing du causes Sd u to decrease sharply, resulting in an increase in the allowed size of the 
couplings cf as seen in Fig. m Approaching d u — 2 causes a rapid increase in Sd u and thus the allowed values of cf 
are again heavily constrained? For du > 2 the quantity Sd u becomes increasingly small, resulting in no bound on the 
couplings cf . 
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FIG. 4: Plot of the allowed cf-c^ parameter space for Bd mixing(left) and B s mixing(right) in the case of two real couplings. 
Black points show regions which agree with the measurement of AMd.s while grey points show additional agreement with the 
measurement of the CP phases <j>d,s- 



We can see from Fig. 2 that there are windows in the allowed parameter space of du — cf . From Eq. (44) we can 



determine r q in terms o\a q as, 

r q = - cos o q ± \jp\- sin 2 a q (61) 

and for p q < 1 the phase a q is constrained to be in the region, 

— 7r — arcsin p q < a q < —n + arcsin p q (62) 

Taking the minimum allowed value of p q we can determine these windows as corresponding to, 

-7T- arcsin p™" + <^ M < arg(^„ ) < -tt + arcsin p" q lm + ^ M (63) 

with $ M = -2j3 s = -0.0409 and = 20 = 0.781, determines these windows to be, 

B d : 1.56 < d u < 1.92 (64) 
B s : 1.80 < d u < 2.00 (65) 

To the inside of these windows, with smaller values of the coupling cf , the unparticle contribution is small enough to 
satisfy the mixing constraint. To the outside of these windows, with larger values of cf , the unparticle contribution 
is larger than the SM contribution. When the phase a q is in the range (— 7r, — ~) the quantity coscr 9 in Eq. (44 1 is 



negative, allowing larger values of r q . This solution corresponds to the case when the unparticle contribution carries 
the opposite sign to the SM and turns over the sign of M\i- 

The black points of Fig. [2] correspond to the region of parameter space allowed by the measurement of the B mixing 
parameters AMj s, while grey points also satisfy the constraint from the measurement of the CP phases 4>d.s respec- 
tively. The left panel of Fig. [2] shows that the allowed region is strongly constrained by the measurement of <f>d in 
addition to the constraint of AM^. This additional constraint restricts the parameter space of c s £ bd considerably. For 
example, for du = § and including the CP constraint we have the improved bound, 

|c£' bd | < 0.024 (66) 




D.2 
S.bd 




FIG. 5: Variation of allowed parameter space of the real coupling cf = cf with scaling dimension du for Bd mixing (left) and 
B„ mixing (right). Black plotted points agree with the CP conserving mixing quantities AMd, s , while grey points also agree 
with the CP phases (j>d,s- 



which is almost three times smaller than the bound from AMg alone. As shown in Eq. (49 1, in the B s system there 
is still a two fold ambiguity in the measurement of the CP phase <p s . This two fold ambiguity is then seen in the right 
panel of Fig. ^ as two distinct pairs of grey regions. The pair of grey regions with larger cf ' hs , and consequently larger 
r a , corresponds to 4>™ p = [-156.90, -106.40]°. This region is clearly disfavoured as it corresponds to the coupling 
cf' bs being outside of the perturbative region. On the other hand the pair of grey regions with smaller cf' bs and r s , 
corresponds to <^ p = [—60.9, —18.58]°. In this case the coupling cf is constrained to be in a narrow band away from 
zero. For example with du — | we have a two-sided bound, 

0.106 < \c s L M \ < 0.23 (67) 

The small cf' bs region also indicates that a value of the scaling dimension in the range 1.22 < du < 1.87 is preferred. 

Despite the two fold ambiguity in <fi s it is clear that values of the scaling dimension in the range du < 2 are preferred 
by the measurement of <p s . Combining these <fi s preferred values for du with the constraints from tpd implies a general 
bound on c L ' bd as, 

\c s L M \ < 0.13 (68) 



2. Two real coupling: cf 7^ 0, cf 7^ 

In this second case we shall allow the left and right unparticle couplings cf and cf to both be real and nonzero. For 
this analysis we shall fix the scaling dimension du = §• The allowed cf — parameter space is plotted in Fig. [ijwith 
black points constrained by AMd. s and grey points further constrained by the CP phase <f>d,s- The allowed region 
extends out along two lines where there is a cancellation between cf and cf . As a result no direct bound can be set 
on the unparticle couplings. 

Starting from the definition of R q = /M^ M and taking the limit m q — ► we arrive at an approximate 

relation, 
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09 „S,bq S.bq 



where the small quantity e q is denned as, 



\M?f M \ 



M B q f B J S du 



(69) 



(70) 



Here it is clear that a large cancellation between the couplings cf and cf. is possible. From Eq. (61 1 we find that 
(p q — l) 2 < r 2 < (p q + l) 2 which leads to, 

- 0.1 x 9 <r q e q < x q (71) 
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with x s — 1 and Xd = A . The solution to these equations is then a parabola in the c£ — c R parameter space described 
by, 



S,bq 



16 
~6 



6e 



If we set e,«0, then these allowed regions follow along the lines, 



S,bq _ , S,bq S,bq 



(72) 



(73) 



in good agreement with the results of the full calculation shown in Fig. [4] 
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FIG. 6: Plot of the allowed |cf |-</>£ parameter space for _B<2 mixing (left) and B s mixing (right) for the case of a single complex 
coupling cf 7^ and = 0. Black plotted points agree with the CP conserving mixing quantities AM^.s, while grey points 
also agree with 4>d,s- 



In this case it is possible to have very large values of the couplings, as long as they lie approximately on the lines 
shown in Eq. (73 1. These solutions with large cancellations are clearly highly fine-tuned. If we consider the case of 



c l = C R then we have no such fine-tuning and we can extract the bounds, 



S,bd 



S,bd 
Z R 
S,bs 



c~£" < 2.5 x 1(T 2 



S,bs o.us , r\ i A 

c L =c R < 0.14 



(74) 
(75) 



from the constraints of pd, s - From the grey regions of Fig. [4] we again see that the measurement of the CP phases 
t>d, s further constrains the allowed parameter space. In the Bd system the allowed region is reduced by more than a 



half. For the case c£ = c R the CP constraint provides a much improved bound, 



S,bd 



S,bd 



< 6.7 x 10" 



(76) 



In the B s system the new CP measurement of (j) s prefers points away from the origin. It is also interesting to notice 
that the case, c S £ hs = c R bs is disfavoured by the latest measurement of 4> s at the 3a level. 

The case of equal left and right couplings is shown in more detail in Fig. [5] where the allowed parameter space is 
plotted as a function of the scaling dimension du- Fig. |5] shows the allowed parameter space for both Bd (left panel) 
and B s (right panel). From the left panel we can see that the allowed parameter space is rather similar to that shown 
in Fig. [2] This time the windows in the plots have disappeared for the case of Bd and have moved in the case of B s 



When we have 



ip we get a new contribution from Cf not present in the case of Sec. 



IIIA1 



ith 



.5 



(I. The 



introduction of Cf changes the overall sign of M^ p and so in this case these holes obey a modified relation similar 



to Eq. (|63| except there is no factor of — tt. This means that for the case of Bd there is no such window and for B s 

as shown in Fig. [5] In this case the impact of the CP phase measurement </><j 
The impact of the recently measured <j) s is much more profound, as 



the window is now at 2 < du < 2.23, 
is to again reduce the allowed parameter space, 
can be seen from the right panel of Fig. [5] where there are no grey points and only black regions. The lack of any 
grey regions in this plot indicates that this case cannot satisfy the 3cr measurement of <p s for any choice of the scaling 
dimension dy. 
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FIG. 7: Plot of the allowed |cf |-|cjj| parameter space for Bd mixing(left) and B s mixing(right) for the case of two complex 
couplings cf 7^ and c% 7^ 0. Black plotted points agree with the CP conserving quantities AMd, s , while grey points also agree 
with (j) d>s . 



3. One complex coupling: cf 7^ 0, cjj = 



In this case we have 



z o = and one complex coupling cf , with the scaling dimension also fixed at du — §. In this 



case the Wilson coefficients again take the simplified form of Eq. (156 h , this time with a complex c S £ bq . As a result the 
phase of A/]^ p is not only determined by the scaling dimension, but also by the phase of the coupling c^' bq . 

Due to the hierarchy nib S> m q the dominant contribution here comes from the Wilson coefficient Of . For the 
purpose of simplicity we can make an approximation of the unparticle contribution as, 



a f?- NP 



g,SM 
12 



where, (5 « ^m 2 



^r.M B jB q B iSd u /M^, with \S\ w 9.45. This then leads to, 

2 



r„ 



arg 



S.bq 



q.SM 
2 



(77) 

(78) 
(79) 



For du = |, we have a q 



2^ 



b s q M , where <f) S L ' bq 



rg ^c^' 69 ^. From Eq. (44) we then have, 



Pq = \l + 2\5\ 



^S,bq 



cos(20 



S.bq 
L 



- z) + \s\ 2 



(80) 



The solutions to this equation are shown in Fig. [6] plotted in the |cf | — </>£ plane. Again we can see from Eq. (62 1 that 
there are holes in the parameter space corresponding to, —38.63° < tj^if" 1 < —6.62° and —65.33 < <j>^ bs < —27.01, 
as determined by Eq. (62 1. Using only the black points of Fig. [6] we can extract a bound on the magnitude of the 
couplings as follows, 



S,bd\ 



< 7.1 x 10" 



|c?' bs | < 0.31 



(81) 
(82) 



From the grey regions plotted in Fig. [6] we can again see that the inclusion of the CP phase constraints leads to a 
much reduced parameter space and improved bounds on the unparticle couplings, 



\c S L M \ < 4.3 x 10- 2 



0.10 < \c S L ' bs 



< 0.31 



(83) 
(84) 



12 





0.2 



a 

-90 -45 45 90 -90 -45 45 90 

i S , bd i E , bd r , - i S.bs i S.bs r i n 

p L = tr L degrees J Q L = Pr LdegireesJ 

FIG. 8: Plot of the allowed |cf |-</>f parameter space for Bd mixing(left) and B s mixing(right) for the case of one complex 
coupling cf = cfi. Black plotted points agree with the CP conserving mixing quantities AMd,s, while grey points also agree 
with cj) d>s . 

In addition we may also constrain the phase of the unparticle coupling c^ bs in the case of B s mixing as, 

-41.77° < <p S £ bs < 25.55° (85) 

4- Two complex couplings: cf 7^ 0, eft ^ 

In this final scalar unparticle case we take two complex couplings, c S £ bq and c^ bq , with the scale dimension fix at 
du = §. Fig. jldisplays the allowed parameter space for this scenario. In a similar way to the case with both couplings 

real, see Fig.Tl the allowed regions extend along the lines \c S £ bq \ ss 5|c^' b9 | and |cf' b9 | « | | c R b9 |- No general bound 
can be placedon the size of the unparticle couplings as large cancellations may occur even for very large values of 
the couplings. These solutions with large cancellations are however very highly tuned and unattractive. Taking the 
subset of this case with the left and right-handed couplings equal, we get the bound, 

\c s £ bd \ = |4' bd | < 0.039 (86) 

\c S L bs \ = \4 bs \ < 0.19 (87) 

In this case the additional constraints from the CP phase measurements of (pd tS have very little effect on the allowed 
parameter space. From the left panel of Fig. [7] it is clear from the clustering of grey points near the origin that the 
CP phase tpd prefer smaller values of |c^' bd |. On the other hand from the right panel we see that the opposite is true 
for |c^' bs | where there are only black points near the origin. 

In Fig. plwe study in more detail the case of equal and complex couplings cf = c^. Fig. |8]shows the experimentally 
allowed |cf]-0f parameter space plotted for scaling dimension du = \- 

The left and right panels of Fig. [8] display exactly the features just described. The grey allowed regions in the left 
panel corresponding to points satisfying both A and <f>d constraints prefer lower values of the coupling in the range, 

| c S,6d| = | c S,M| < Q Q24 

From the right panel we can see that the phase <p s indeed prefers regions with a larger coupling in the region, 

0.06 < |cf' bs | = \c s / s \ < 0.18 (89) 
In Fig. 0we again see holes in the |cf | — 4>f parameter space. In the case of cf = c^ we h ave a a = 2<\> S L bq 



a,rg(Sd u Y~ 4>q M , where the extra ir is from a change in sign of A/^ p relative to the case in Sec. Ill A3 Using Eq. (62 1 
we find that these holes correspond to, 

24.67° < <j) S L M < 62.99° (90) 

51.37° < <j) s £ bs < 83.38° (91) 
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FIG. 9: Plot of the variation of the allowed c\ parameter space with the scaling dimension du for Bd mixing(left) and B s 
mixing(right) in the case of a single real coupling c\ 7^ and c\ = 0. Black plotted points agree with the CP conserving 
mixing quantities AMd, s , while grey points also agree with <f>d, s - 



B. Vector unparticles 



Now let us consider the case of vector unparticles. In general the contributions from vector unparticles to meson 
mixing is larger than from scalar unparticles due to the enhancement of the Wilson coefficients by the factor A^/Mg . 
As a result the experimentally allowed parameter space for vector unparticles will in general also be suppressed by 
this same factor. 



1. One real coupling: 7^ 0, Cr = 



First we shall consider the case of a single real vector unparticle coupling, c V rf hq . Assuming c^ bq — we then choose 
to allow the scaling dimension du to vary and study the experimentally allowed parameter space. This case is rather 
similar to that of scalar unparticle case discussed in Sec. |III A 1| The allowed regions of parameter space for this case 
are shown in Fig. [9]with black points indicating the regions of parameter space allowed by the measurement of AM<j )S , 
while grey points indicate in addition agreement with measurements of the CP phase <pd,s- Fig. [9] again shows that 
the measurement of the CP phases (fid.s improves greatly on the constraints set by AMj s . For example with du = § 
the constraint from AMd, s alone gives, 

\c V rf bd \ < 6.8 x 1(T 4 (92) 
\c V rf bs \ < 2.2 x 1(T 3 (93) 
Adding the restriction from the measurement of the CP phases cj>d s we get the improved bounds, 

\c V rf bd \ < 2.4 x 1(T 4 (94) 
1 x 1CT 3 < \c V rf bs \ < 2.2 x 1CT 3 (95) 

Looking at the grey allowed regions of the right panel of Fig. [9] it is clear that the recently measured phase in B s 
mixing, <p s , prefers values of the scaling dimension in the range 1.22 < du < 1.96. With the scaling dimension du in 
the range 1 < du < 2 we may extract the general bounds, 

\cl M \ < 1.3 x irr 3 (96) 

2.4 x 1(T 4 < \cl' bs \ <3.1 x 1(T 2 (97) 



Two real couplings: c£ 7^ 0, c\ 7^ 



In this case we fix the scaling dimension du = § with two real nonzero couplings c£ bq and c^ bq . 

V,bq V.bq 



Fig. 



10 



shows 



the allowed regions of the c L ' -erf parameter space. Here the unparticle contribution to the AF = 2 effective 
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FIG. 10: Plot of the allowed c^-Cb, parameter space for Bd mixing(left) and B s mixing (right) in the case of two real couplings 
c\ 7^ and c)j 7^ 0. Black plotted points agree with the CP conserving mixing quantities AMd, s , while grey points also agree 
with (j) d>s . 



Hamiltonian satisfies the approximate relation, 



{c VM f + {c VM f 



ToS,bqS,bq 
iz C L C R 



ML 



1 c 9 A 2 
1V U 



(98) 



where the small quantity e q is as defined in section III A 2 Again solutions extend along two lines of the approximate 
form. 



-L 

n V,bq 



-0.3 c 
-3.3 c 



V.bq 
R 

V,bq 
R 



(99) 

(100) 



as can be seen in Fig. |10| In general no bound can be placed on the size of the unparticle coupling in this case due 
to the possibility of large cancellations. Let us then again consider the simplified case with left and right unparticle 
couplings set equal to each other, c^ hq — c^ bq . First using the constraints from AMd, s we may set the bounds, 

I V,bd 



\cl M \ < 2.7 x 10- 4 



V.bs 



'R 

V,bs 
"R 



< 9.2 x 10~ 



(101) 
(102) 



The grey plotted points of Fig. 10 show the effect of the additional constraint from the CP phase measurements of 
(j>d, s - The allowed parameter space is once again much reduced by the inclusion of these CP constraints. The CP 
phase (j>d improves the above bound on c^ bd to, 



V,bd, 



\cl' bd \ < 9.6 x 10- 5 . 



while the CP phase constraint of </> s improves the above bound on c^ hs to 



(103) 



when du = § • 



4.2 x 10- 4 < \cl M \ 



< 9.2 x 10" 



(104) 



The case of 



1 y 



is displayed in Fig. 



11 



niBi 



which very much resembles the case of section 
0. Here with two identical real couplings = the allowed parameter space is less than 



with just one real 



coupling ^ and c v R 

half that found in section |III B 1| As such all bounds on the couplings found in section |III A 1| may be restated for 
this case by simply multiplying by 0.42 for B s and 0.4 for Bd- 
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FIG. 11: Plot of the variation of the allowed parameter space with the scaling dimension du for Bd mixing(left) and B a 
mixing(right) in the case of one real coupling c\ = Cr- Black plotted points agree with the CP conserving mixing quantities 
AMd,s, while grey points also agree with 4>d,s- 




-45 
<|) L ' [degrees] 




FIG. 12: Plot of the allowed \c\\~4>l parameter space for Bd mixing(left) and B s mixing(right) in the case of one complex 
coupling c[ 0, c\ — 0. Black plotted points agree with the CP conserving mixing quantities AMd, s , while grey points also 
agree with <f>d,s- 



3. One complex coupling: 



Next we consider the case of a single complex coupling and fix the scaling dimension to be du 



This case is 



similar to that analyzed in section III A3 The allowed parameter space for this case is shown in Fig. 12 The black 
points in Fig. 12 show the regions of parameter space which agree with experimental measurements of AM^ S while 
grey points indicate the additional agreement with the CP phases <pd,s- From the black points we can extract the 
bound, 



S,bd, 



S,bs, 



< 7.2 X 1CT 4 

< 2.9 x 1CT 3 



(105) 
(106) 



It is clear from the grey plotted points of Fig. 12 that the CP phase constraints have an important role to play in 
this scenario also. From these grey regions we may again extract improved bounds in the magnitude of the couplings, 



9.7 x 10 -4 < 



S,bd, 



S,bs, 



< 4.4 X 10~ 

< 2.9 x 10~ 



(107) 
(108) 



1G 



In addition we may also constrain the phase of the unparticle coupling c^ bs as 



41.80° < 4>l M < 25.52° 



(109) 
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FIG. 13: Plot of the allowed |cj*|-|c)j| parameter space for Bd mixing (left) and B a mixing (right) for the case of two complex 
couplings. Black plotted points agree with the CP conserving mixing quantities AMd, a , while grey points also agree with (f>d,s- 



4- Two Complex Couplings: c£ 7^ 0, c)j ^ 



In this final scenario we take the general case of two independent complex couplings c£ bq and c^ bq , and fix dy 



3 

The experimentally allowed parameter space for q = d,s are shown in Fig. [13] again showing the allowed regions from 
both AMd s and (f>d s constraints. This case is closely related to that of section III A 4 where we found that the allowed 



regions followed along two lines. Again here the allowed regions approximately follow the lines |c^ 09 | = 0.3 |c^ &9 | and 



3.3 |c^' b9 | and as a result no general bound can be set in this case either. From the grey regions plotted in 



Fig. 13 we can see that the CP phases <pd yS prefer smaller |c y ' M | and larger |c 



V.bsl 



Let us consider the case when the two couplings are equal. This case is plotted in Fig. [14J in the plane of \c\ 
The CP conserving Bd, s mixing constraints shown by the black plotted points may provide the bounds, 



S,bd, 



S,bs, 



< 2.9 X 10" 

< 1.2 x 10" 



(110) 
(111) 



When we include the CP phase constraints of <f>d, s we get the improved limits, 



|cf ,M | < 1.8 x 10- 4 
4xl0- 4 < \c S L ' bs \ <1.2xl0- 3 



(112) 
(113) 



In this case we may again constrain the phase of the unparticle coupling c^ bs as 



41.80° < (j) V £ bs < 25.52° 



(114) 



IV. CONCLUSIONS 



In this work we have studied the contribution of the recently suggested unparticle physics to B meson mixing. Our 
main aim was to study the impact of the CP phase associated with Bd, s mixing. In particular we have been interested 
in the impact of the recent 3cr evidence for new physics found in the phase of B s mixing. The phase of Bd mixing is 
rather close to the standard model expectation and as such these two phases provide very different constraints on the 
allowed phase of new physics. 
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FIG. 14: Plot of the allowed \c£\-(pL, parameter space for B d mixing(left) and B s mixing(right) for the case of one complex 
coupling cl = c\. Black plotted points agree with the CP conserving mixing quantities AMd iS) while grey points also agree 
with cj) d>s . 

For our numerical analysis we have chosen to study the cases of scalar unparticles and vector unparticles separately. 
Each case has further been subdivided into a number of scenarios depending on the type unparticle couplings. In 
each scenario we have been careful to understand and discuss the specific characteristics of the allowed parameter 
space and, in particular, to understand the impact of the CP phase constraints on this parameter space and, where 
possible, to set bounds on the unparticle couplings. Our main goal here has been to identify the impact of the phase 
constraints on the allowed parameter space in each scenario and, in particular, that of the recently measured phase 
of B s mixing, <fi s . A summary of the bounds placed on the unparticle couplings from the mixing parameters AA/^s 
alone and for AMj s together with are shown in Table ITJ 



Scalar unparticles 


AM d 


A M s 


AM d & 4> d 


A M s & (j) s 


cf G R, eg = 


\cf bd \ < 0.067 


\c b L ' ba \ < 0.23 


\c s L M \ < 0.024 




0.106 < \c b L ' b3 \ < 0.23 


cf, eg G R 


no constraint 


no constraint 


no constraint 




no constraint 


eg = eg G R 
eg G C, eg = 


|cg M j < 0.025 
cg M < 0.071 


\cf ba \ < 0.14 
|cg' 6a i < 0.31 


|cg' M | < 0.0067 
|cg' M | < 0.043 




Excluded 
0.10 < |cf' 6s | < 0.31 


eg, eg G C 


no constraint 


no constraint 


no constraint 




no constraint 


eg = eg G C 


|cg M j < 0.039 


|c£' 6s | < 0.19 


|cg' M | < 0.024 




0.06 < \c s L ' bs \ < 0.18 


Vector unparticles 




AM S 


AM d & 4> d 


AM S & 4> s 


cl G R, c v R = 


|c£ M | < 6.8 x 10" 4 


\c*' bs \ < 2.2 x 10" 3 


\c^ M \ < 2.4 x lCT 


4 


1 x 10 -3 < |c^ 5s | < 2.2 x 10" 3 


cl, clem 


no constraint 


no constraint 


no constraint 




no constraint 


cl = cleR 


|cg' bd | < 2.7 x 10" 4 


\c s L ' ba \ < 9.2 x 10- 4 


|cg' M | < 9.6 x 10" 


5 


4.2 x 10~ 4 < |cg' 6s | < 9.2 x 10~ 4 


cl G C, cl = 


cg M | < 7.2 x 10" 4 


cg'H < 2.9 x 10" 3 


c S,6d < 4A x 1Q - 


4 


9.7 x 10" 4 < \c s L ' bs \ < 2.9 x 10" 3 


cl, clec 


no constraint 


no constraint 


no constraint 




no constraint 


cl = cl G C 


|cg M | < 2.9 x 10" 4 


|c£' 6s | < 1.2 x 10" 3 


|cg' M | < 1.8 x 10" 


4 


4 x 10 -4 < |cg bs | < 1.2 x 10~ 3 



TABLE I: Summary of bounds on unparticle couplings with scaling dimension du = § from 95% CL. measurement of AMd,s 
alone and also from a combination of the 95% CL. measurements of both AM d<s and <f)d,s- 

The main conclusion of this work is that the CP phase constraints do indeed improve markedly on the CP conserving 
bounds set on unparticle couplings in every scenario studied by a factor of 2 ~ 4, see table |T] Of particular interest is 
the case of scalar unparticles with real and equal couplings c S £ bs = c^ bs for the B s system. Using only the measurement 
of AM S we have a bound in the region of 0.14, but by adding the constraint from the 95% C.L. measurement of <j) s 
this case is excluded. 
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It should be noted that the combination of D0 and CDF results for (j> s were made without full knowledge of the likelihoods, 
which could have an effect on the significance of this SM deviation. These are now being made available and a new 
experimental combination is due to be released soon. Until that time this present hint of new physics is the best available 
and remains a very exciting prospect. 



